Decoherence of a particle in a ring 
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Abstract. - We consider a particle coupled to a dissipative environment and derive a perturbative 
formula for the dephasing rate based on the purity of the reduced probability matrix. We apply 
this formula to the problem of a particle on a ring, that interacts with a dirty metal environment. 
At low but finite temperatures we find a dephasing rate oc T 3 / 2 i and identify dephasing lengths 
for large and for small rings. These findings shed light on recent Monte Carlo data regarding the 
effective mass of the particle. At zero temperature we find that spatial fluctuations suppress the 
possibility of having a power law decay of coherence. 
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Introduction. — The problem of dephasing of a par- 
ticle coupled to a dissipative environment at tempera- 
ture T, and in particular in the limit T — > has fas- 
cinated the mcsoscopic community during the last two 
decades [1-7]. It has been shown [8,9] that the Caldeira- 
eggett (CL) framework [10, 11] can be generalized and 
that the proper way to characterize the environment is 
by its form factor S(q,u>). Application of the Feynman- 
Vernon formalism [8, 9] and a semiclassical analysis have 
shown that an interference amplitude P v decays with time 
as P v = exp(— Pip(t)) with 



p v (t) 



t f f S(q,u;)P(-q, 



-u) 



(1) 



where the integration measures over the wavevector and 
the frequency are d 3 g/(27r) 3 and dui/2n respectively. The 
interference suppression factor P v is known in the litera- 
ture as the dephasing [2] or as the decoherence factor [25], 
and in the present work we show that it reflects loss of 
purity. In the semiclassical treatment S(q, ui) is the sym- 
metrized form factor of the environment and P(q, ui) is 
the classical symmetric power spectrum of the motion. It 
has been conjectured using ad-hoc argumentation [12] (see 
also [2,13]) that "the correct" procedure is to use non- 
symmetrized spectra. One of our aims is to provide a 
proper derivation for a corrected Eq.(l). 

During the last decade the study of a particle in a 
ring coupled to a variety of environments, has become a 
paradigm for the study of ground state anomalies [14-18]. 
Besides being a prototype model problem it may be re- 
alized as a mcsoscopic electronic device, and it is also of 
relevance to experiments with cold atoms or ions that are 



trapped above an "atom chip" device [19-21], where noise 
is induced by nearby metal surfaces. A significant progress 
has been achieved in analyzing the equilibrium properties 
of this prototype system, in particular the dependence 
of the ground state energy on the Aharonov Bohm flux 
through the ring. 

In the present work we define "dephasing" as the pro- 
gressive loss of purity and find a consistent revised form of 
Eq.(l) that is valid beyond the semiclassical context. We 
apply this result to the model of a particle on a clean ring 
that interacts with a dirty metal environment. At finite 
temperature we identify the dephasing rate T v = p v (t)/t, 
that vanishes at zero temperature. At T=0 we find that 
only in the CL-like limit of our model there is still slow pro- 
gressive spreading {p v (t) ~ hit) which suggests a power 
law decay of coherence. Our results shed new light on re- 
cent Monte Carlo data for the temperature dependence of 
mass-renormalization [22]. 

Purity. — Our starting point is the most natural def- 
inition for the dephasing factor as related to the purity 
trace(p 2 ) of the reduced probability matrix. The notion 
of purity is very old, but in recent years it has become 
very popular due to the interest in quantum computa- 
tion [23] . Assume that the state of the system including 
the environment is ^ pn , where p and n label the basis 
states of the particle and the bath respectively. Tracing 
the environment states n defines a reduced probability ma- 
trix [p aya ]p,p> = J2 n ^pn^p'n and tne purity is then mea- 

Assum- 



sured by the dephasing factor P v = Wtrace^^J 
ing a factorized initial preparation as in the conventional 
Feynman- Vernon formalism, we propose the loss of purity 
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(P v < 1) as a measure for decoherence. A standard reser- 
vation applies: initial transients during which the system 
gets "dressed" by the environment should be ignored as 
these reflect renormalizations due to the interactions with 
the high frequency modes. Other choices of initial state 
might involve different transients, while the later slow ap- 
proach to equilibrium should be independent of these tran- 
sients. In any case the reasoning here is not much different 
from the usual ideology of the Fermi golden rule, which is 
used with similar restrictions to calculate transition rates 
between levels. 

Consider then a factorized initial preparation 
*<°) 

^ pn 



Sn,: 



yp,p <Jn,n , 



so that within perturbation theory 



all fy pn arc small except for ^ Po n - We can relate P v 
to the probabilities Pt(p, n\p , n ) = l^pn] 2 to have a 
transition from the state Ipoj^o) to the state \p,n) after 
time t. To leading order we find 

P v = ft(Po,«obo,«o) (2) 
+ E Pt(p,n \po,n ) + ^2 Pt(Po,n\p 0> n ) 



P^Pa 



The first term in Eq.(2) is just the survival probability 
Psurvivai of the preparation. The importance of the two 
other terms can be demonstrated using simple examples: 
For an environment that consists of static scatterers we 



have P., 



< 1 but P w 



1 thanks to the second term. 



For a particle in a ring that interacts with a q = environ- 



mental mode P., 



< 1 but Pu> = 1 thanks to the third 



term. Using n Pt{p, n\po, uq) = 1 we finally obtain 



P<p 



1-P V = 



E 



pjtpo nutria 



Pt{p,n\p 0l n ) 



(3) 



This result has the form of a Fermi-golden-rule (FGR), 
i.e. it is the probability that both the system and the 
bath make a transition. This differs from the usual FGR 
treatment [2] in which terms like Pt(po,n ^ no\po,no) 
are included. In the problem that we consider in this 
paper we can calculate P v using a dqdoj integral as in 
Eq.(l). In many examples the w = transitions have 
zero measure and therefore P v is practically the same as 
Pt{po, «o|Po, n o)- Otherwise one has to be careful in elim- 
inating those transitions that do not contribute to the de- 
phasing process. Anticipating the application of Eq.(l) 
for the calculation of the dephasing for a particle in a 
ring, the integration over q becomes a discrete summation 
where the <j=0 related component should be excluded. It 
is implicit in the derivation of Eq.(l) from Eq.(3) that at 
the last step a thermal average is taken over both Uq and 
POj though in general one may consider non-equilibrium 
preparations as well. 

Dephasing formula. — We would like to apply our 
revised FGR Eq. (3) to the general problem of a particle at 
position R coupled to an environment with electronic den- 
sity n(r, t). It is implicit that the particle also experiences 
an external potential that defines the confining geometry. 



A Hamiltonian Tio of the particle in the confined geome- 
try defines the states and eigenstates via Tto\p) = E p \p). 
For definitcness we use the Coulomb interaction, though 
any other interaction may be used, hence the particle - 
environment interaction is 



(4) 



K.t = J d\p(v) J dV^^ = J d\p{v)U{v) 

where p(r) = 6(r — R(t)), with R(t) the position opera- 
tor of the particle in the Hciscnbcrg (interaction) picture. 
Our FGR with P t (p, n\p , n ) = |(p, n\ f Q V int dt'\p , n )\ 2 
(using h = 1 units) yields 



Pr 



2 E E 

p(¥=Po) n(=£n a ) 
Jl +H 



t <-t 

dt' / dt" 
Jo 



(p \p(v" 1 t")\p)(p\ P (r\t')\p ) 
(n \U(r",t")\n}{n\U(r',t')\n ) (5) 

The double time integral can be written as a dqdu integral 
over Fourier components. For this purpose we define the 
form factor of the fluctuations (as seen by the particle) : 



S(q,u)= / d 3 r / dT(U(r',t')U(r,t))e luJT - lcl - 



(6) 



with thermal average replacing the no state expectation 
value. S(q, ui) is related to the dielectric function of the 
environment s(q,uj) via the fluctuation dissipation theo- 
rem 
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S(q, uj) = — 5-Im 



-1 



1 



=-w/T 



(7) 



_e{q,u}) 

In the scmiclassical formulation one replaces the operator 
R(t) by the classical trajectory R c i(t), and consequently 
in Eq.(5) the particle-related part of the integrand is re- 
placed by a classical two-point correlation function of the 
type (f(R c i(t"))f(R c i(t'))}. In the quantum context the 
particle dependent part of Eq.(5), after Fourier transform, 
leads to the following definition for the power spectrum of 
the motion: 



Also here, close to equilibrium condition, a thermal aver- 
age should replace the po state expectation value. It is 
important to realize that this definition, as well as Eq.(7), 
imply that non-symmetrized spectral functions should be 
used. Our main interest is in very low temperatures, so 
we set for presentation purpose po = 0. Then we get 



P(g, W )=X>|e 

p#0 



iq-Rl 



S(u - E p ) 



(9) 



Using the above definitions for S(q, uS) and P{q 1 u>) we can 
re- write 



t ft 



p v = I dt' I dt" I // S(q,u/)P(-q,-u") 

10 JO JqJJu',ui" 

xe -i(u'- U ")(t"-t') _ ( 1Q ) 
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For practical calculations or for aesthetic reasons we prefer 
to use soft rather than sharp cutoff for the time integra- 
tion. Then we get 



with 



t 



S(q,u)P(-q, 



(2/t) 



(1/i) 2 + (lo-lo') 2 



This result can be cast into the form of Eq.(l) provided 
P(q,u>) is re-defined as the convolution of Eq.(8) with 
the kernel in the square brackets, which is like time- 
uncertainty broadened delta function. Eq.(ll) is our re- 
vised form of Eq.(l); it provides the dephasing factor 
P v = exp(—p v ) for a general particle-environment inter- 
action. 

Dephasing rate. — At finite temperatures, if t is 
larger compared with dynamically relevant time scales, 
and in particular t > 1 /T, we can replace the square 
brackets in Eq.(ll) by 2ir5(ui— ui'). Consequently we get 
linear growth p v 



Y v t with the rate 



S(q,u)P(-q,-u) 



(12) 



Following standard arguments one conjectures that the 
long time decay is exponential, i.e. P(t) = exp(— T^t), as 
in the analysis of Wigncr's decay. In terms of the dielectric 
function e{q,w) we obtain the following general result: 



r. 



E 

p#0 
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-Im 



i 



e(q,-E p ) 



2\{p\c 



iq-Rl 



i E p /T_ 1 



(13) 



Our assumption t^>(l/T) implies that (13) can be 
trusted only if T v <C T which implies a weak coupling con- 
dition (see below). 

Dirty metal. — So far we kept the derivation general, 
without specifying either the particle states \p) or the di- 
electric function e(q,uj). We consider now a particle of 
mass M on a ring of radius R so that Ho = — (2MR 2 )~ 1 d 2 
where 9 is the angle variable. The particle eigenstates 
are then \p m ) oc c" n6 with energy eigenvalues E m = 
m 2 /(2MR 2 ). We study the effect of low frequency fluctu- 
ations (\q\ < l/£, \lu\ < uj c ) due to a dirty metal environ- 
ment for which e(q,ui) = 1 + 47r<r(— iuj + Dq 2 )" 1 , where a 
is the conductivity, D is the diffusion constant, and £ is the 
mean free path. Below we identify the renormalized value 
of the high frequency cutoff u> c as the classical damping 
rate j r = 2ira/M£ 2 , where the dimensionless interaction 
strength is a = e 2 /(87r 2 a£) = 3/(8(k F £) 2 ) and k F is the 
Fermi wavevector. We first consider the case of a large ring 
with r = R/l ^> 1. Using the Fourier expansion [16, 17] 



£ /" e -iq-(R(0)-R(0'))fZ[ 
Jq <f 



1 



4 r 2 sin 2(^_91) + 1 



m{6-0') 



(14) 



— In ■ 

izr 

o, 



1 < m < r 
otherwise 



(11) we have 



f ^ K0 | e -^ R |p m )| 2 = ia H 

Jo d ^ 



(15) 



(16) 



and therefore 



T v = 2ira ^ < 



i ^e-JtZ]- ~ 2lTaT E am ( 17 ) 

0<|m|<r e f£ 



where r off ee min{r, {2MR 2 T) 1 / 2 , {2MR 2 oj c ) 1 / 2 } is deter- 
mined by the conditions m < r and E m < T and E m < uj c . 
At high temperatures, T > uj c , r eff is temperature inde- 
pendent and therefore T v oc T, while at low temperatures 
[but still T > 1/(2M_R 2 )] we get, as shown schematically 
in Fig. 1, 



T v = AaTV2M£ 2 T\ \n\ / 2Ml 2 T\ ~ T 3/2 | lnT| (18) 

From these results it follows that the self consistency re- 
quirement r v -c T, as discussed after Eq.(13), is globally 
satisfied for any temperature if a <C 1, or cquivalently if 
kp£ 3> 1; in the regime of Eq.(18) the constraint T v <C T 
is satisfied also with stronger a. We note that if the m=0 
Fourier component were included in the summation, then 
the low T form would change to T v oc T, in contrast with 
the proper result Eq.(18). 

Though the derivation of Eq.(ll) refers to a "one par- 
ticle" , it turns out that the treatment of dephasing in 
the many body problem is not much different. Namely, 
the effect of the " Pauli principle" in the low temperature 
Fermi sea occupation can be incorporated via an appro- 
priate modification of the cutoff scheme for the spectral 
functions involved. On the heuristic level it involves a T 
dependent momentum cutoff as in [12], while more re- 
cently it has been formulated and established using more 
advanced methods [24]. 

Zero temperature. — We would like to discuss the 
"zero temperature" regime. If the temperature were ex- 
tremely low, such that T <C A where A ~ l/(Mi? 2 ) is 
the ground state level spacing, we could treat the prob- 
lem using a "two level approximation" , which is a very 
well studied model [25]. We do not further discuss this 
regime. From here on we assume T A. Thus we can 
treat the dio integration as if the levels of the ring form 
a continuum. But we still can define "zero temperature" 
as such for which the practical interest is in the time in- 
terval t <C which can be extremely long. Then one 
realizes that Eq.(ll) gives a non-zero result even at "zero 
temperature" : 



P v 



a 



E a - 



In- 



E„ 



(iA) 



(19) 
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can be parameterized as 



Fig. 1: Illustration of the dependence of the dephasing rate V 
on the temperature T. The dephasing rate is well defined 
for t > (1/T), and hence the self consistency requirement is 
r <C T. This condition is demonstrated by a comparison 
with the dashed line. The illustration assumes weak coupling 
a < 1 and large rings ar 2 > 1 so that the energy cutoff is 
7,. = 2-KOt/Mi 2 > A = 1/MR 2 . For extremely low tempera- 
tures, such that T is smaller compared with the spacing A, the 
probability to excite the system is exponentially small and the 
familiar two-level modeling becomes applicable. 



where uj c is the high frequency cutoff of the environmental 
modes. Assuming I <C R we get after a transient: 



af 



dq In 



1 



In 



MiOr 



a (20) 



where for clarity of presentation we converted the m sum- 
mation into a dq = (l/R)dm integral. Accordingly we de- 
duce that for r ^> 1 coherence is maintained if kpt 1. 
This should be contrasted with the CL limit (£ — > oo) 
where the integral has a singular q ~ contribution from 
the lowest fluctuating mode (m = 1), and consequently 
p v « 2ar 2 log(w c t), which is a well known expression 
[6,15]. But once r < 1 the condition for using our per- 
turbation result is replaced by ar 2 -C 1, which is an I in- 
dependent condition. In this CL limit the quantization 
of the energy spectrum is important and the renormal- 
ized cutoff frequency becomes tu c ~ A instead of u> c ~ j r . 
Accordingly, in the latter circumstances, the time during 
which the log spreading prevails diminishes. 

Effective mass. — It can be shown [9] that the 
mass rcnormalization in the incrtial (polaronic) sense is 
AM = rj/u) c . However in recent works [16-18] the mass 
rcnormalization concept appears in a new context. The 
free energy F(T, $) of a particle in a ring is calculated, 
where $ is the Aharonov Bohm flux through the ring. 
Then the coherence is characterized by the "curvature", 
which is a measure for the sensitivity to $. The curvature 



o 2 f 



*=0 



M*R 2 



f(M*R 2 T) 



(21) 



where in the absence of environment M* = M is the bare 
mass of the particle, and the T dependence simply reflects 
the Boltzmann distribution of the energy. In the presence 
of coupling to the environment M* > M and M* depends 
on both a and T. At T = 0, for fixed a < 1, Monte 
Carlo data show [22] that the ratio M* /M is independent 
of the radius provided r > r c , where r c is a critical ra- 
dius. As the radius becomes smaller compared with r c , 
the ratio M* /M rapidly approaches unity. In the regime 
of "large R" the mass renormalization effect diminishes 
with the temperature and depends on the scaled variable 
RT, while for "small R" the ratio M*/M depends on the 
scaled variable i? 4 T. The natural question is whether we 
can shed some light on the physics behind this observed 
temperature dependence. Making the conjecture that the 
temperature dependence of M*/M is determined by de- 
phasing it is natural to suggest the following measure of 
coherence: 



2nadMR z T 



(22) 



where a is an average value of a n . Eq.(22) de- 
scribes dephasing at the time t = 1/A ctf , where 
A off ~ r cff x (MR 2 )" 1 is the energy scale that character- 
izes the "effective" transitions; hence the variable x mea- 
sures the level sharpness. For a dirty metal with £ <§; R 
the typical value of the Fourier components is a ~ 1/r as 
implied by Eq.(15). On the other hand for a dirty metal 
with I ^> R there is only one effective mode with a ~ r 2 . 
Accordingly we get the RT and the R 4 T dependence re- 
spectively, in agreement with the Monte Carlo simulations. 
The condition x < 1/2 can serve as a practical definition 
for having coherence. It can be translated either as a con- 
dition on the temperature, or optionally it can be used 
in order to define a coherence length that depends on the 
temperature. The conjecture is that M* /M is a function 
of x. 

Summary. — In this paper we derive a new perturba- 
tive expression for the dephasing factor P v (t) and apply 
it to a particle in a ring coupled to fluctuations of a dirty 
metal environment. We find that the dephasing rate van- 
ishes at T = 0. We also define a coherence criterion that 
identifies a dephasing length. The latter diverges as T _1 
for large radius and as T^ 1 / 4 for small radius, in consis- 
tency with Monte Carlo data on mass renormalization. 
The renormalizcd mass is an equilibrium property which 
affects temporal correlation functions. But we see that it 
reflects nonequilibrium features of the dynamics which are 
expressed in the dephasing factor calculation. We find this 
relation between equilibrium and nonequilibrium scales an 
intriguing phenomena. 
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